Draft version January 24, 2009 

Preprint typeset using I4TgX style emulateapj v. 10/09/06 



GRS 1915+105 IN "SOFT STATE": NATURE OF ACCRETION DISK WIND AND ORIGIN OF X-RAY EMISSION 

YOSHIHIRO UEDA 1 , KAZUTAKA YAMAOKA 2 , RONALD REMILLARD 3 

Draft version January 24, 2009 

ABSTRACT 

We present the results from simultaneous Chandra HETGS and RXTE observations of the microquasar GRS 
1915+105 in its quasi-stable "soft state" (or State A) performed on 2007 August 14, several days after the 
state transition from "hard state" (State C). The X-ray flux increased with spectral hardening around the middle 
of the Chandra observation, after which the 67 Hz quasi periodic oscillation (QPO) became significant. The 
HETGS spectra reveal at least 32 narrow absorption lines from highly ionized ions including Ne, Mg, Si, S, Ar, 
Ca, Cr, Mn, Fe, whose features are the deepest among those ever observed with Chandra from this source. By 
fitting to the absorption line profiles by Voigt functions, we find that the absorber has outflow velocities of ss 150 
and «500 km s" 1 with a line-of-sight velocity dispersion of «70 and w200 km s" 1 for the Si XIV and Fe XXVI 
ions, respectively. The larger velocity and its dispersion in heavier ions indicate that the wind has a non-uniform 
dynamical structure along the line-of-sight. The location of the absorber is estimated at ~ (1 -3) x 10 5 r g (r g 
is the gravitational radius) from the source, consistent with thermally and/or radiation driven winds. By taking 
into account narrow spectral features detected with Chandra, the continuum spectra obtained with RXTE in 
the 3-25 keV band can be well described with a thermal Comptonization with an electron temperature of ~ 4 
keV and an optical depth of w 5 from seed photons from the standard disk extending down to (4 — 7)r g . In 
this interpretation, most of the radiation energy is produced in the Comptonization corona, which completely 
covers the inner part of the disk. A broad (ler width of w 0.2 keV) iron-K emission line and a smeared edge 
feature are detected, which can be explained by reflection from the accretion disk at radii larger than 400r g . 
Our data do not require the presence of an extremely blurred iron-K line or disk emission originating from radii 
smaller than w 4r g in this state. 

Subject headings: accretion, accretion disks — stars: individual (GRS 1915+105) — techniques: spectroscopic 
— X-rays: stars 



1. INTRODUCTION 

Microquasars are key objects to investigate the physics of 
accretion onto a black hole in relation with relativistic jet. 
GRS 1915+105 is the most important, archetype microquasar, 
from which the first super luminal motion in our G alaxy was 
detected (for review, see iFender & Bellon! 12004). It' is an 
X-ray binary consisting of a M-K III star and a black hole 
of 14 ±4 M Q with an inclination angle of 66°-70° at a dis- 
tance of « 12 kpc. The luminosity of GRS 1915+105 reaches 
close to the Eddington limit, indicating critical mass accretion 
takes place. With this reason, GRS 1915+105 exhibits unique 
properties distinct from canonical black hole binaries (BHBs). 
In particular, it sometimes shows dramatic temporal/spectral 
variations ("oscillation") occurring in quasi-regular cycles, 
most probably due t o thermal instability in the accretion disk. 
iBelloni et alJ ((2000) phenomenologically classified the X-ray 
states into three, States A, B, and C. GRS 1915+105 spends 
most of time in State C where the X-ray spectrum is relatively 
hard. The underlying physics of these complex behaviors and 
the relation to canonical BHBs is far from being fully under- 
stood, even though it is a key question to establish a general 
solution of accretion onto black holes at various mass accre- 
tion rates. 

Outflow by a disk wind is one of key ingredients in under- 
standing the dynamics of an accretion flow in dis k systems, in- 
cluding X-ray binaries and active galactic nuclei. Kota ni et"ai1 
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(2000) discovered absorption line features of iron-K ions in 
the X-ray spectra of GRS 1915+105 with ASCA, revealing 
the presence of a significant amount of highly ionized plasma 
in this system. Similar features were reported in many other 
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which are now recognized to be a 
common feature in accreting sources over a wide range of 
luminosity. The origin of such ionized gas has been estab- 
lished to be a disk wind by the precise measurement of the 
outflow velocity with Chandra. The physical mechanism to 
produce the disk wind is not established yet. There are at least 
3 possibilities that are not exclusive one another: thermally- 
driven wind (LBe gelman et al. 1983), radiation driven wind 
(e.g., iProga. Stone. & Kallmanl l2000h . and/or magnetically 
driven wind. Based on a simple argument from an observed 
column density and an ionization parameter, the location of 
the wind has been estimated at 10 4_5 r g (r g = GM/c 2 is the 
gravitational radius where G, M, c is the gravitational con- 
stant, mass of the black hole, and light velocity, respectively.) 
from the source, consistent with thermally and/or radiation 
driven in all the ca ses (including GRO J 1655^-0) except for 
iMiller et"ai1 (120061) . who deduced a much smaller radii from 
Chandra data of GRO J 1655-40, leading to their conclusion 
that a magnetically drive n wind i s the o nly solution. This ar- 
gument is questioned by Netzer (2006), however, making the 
discussion still controversial. 

Detailed studies of the disk wind in GRS 1915+105 is par- 
ticularly important to understand the role and origin of the 
outflow at high mass accretion rates. Among the different 
epochs observed with ASCA, the absorption features were the 
most prominent (equivalent widths of 30-50 eV for H- and 
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He-like Fe ions) in 1994 and 1995 when the spectrum is soft 
and the 1-10 keV flux is relatively faint (—0.3 Crab). We 
recognize that in these ASCA observations the source was in 
quasi-stable "soft state" (State A), which is found to be quite 
rare case in the long term variability from the RXTE/ASM 
over >10 years. By contrast, the iron-K absorption-line fea- 
tures were much weaker and even undetectable when the 1- 
10 keV spectrum was hard, corresponding to State C (or "hard 
state"). This can be explained by the increasing photon flux 
responsibl e for photo-ioni zation which almost fully ionized 
the ions. iLee et al] d2002l) performed a HETGS observation 
of GRS 1915+105 in State C and detected absorption lines 
from highly ionized Fe ions. However, the features are much 
weaker, with equivalent widths of only 10 eV, compared with 
those found in State A with ASCA. Because of the weakness 
of these features there still remain uncertainties in the basic 
properties of the plasma, such as the bulk velocity and its dis- 
persion (or kinetic temperature). 

In this paper, we report the results of high energy resolu- 
tion spectroscopy of GRS 1915+105 in "soft state" (State A) 
performed simultaneously with Chandra HETGS and RXTE 
in 2007 August. The first purpose is to reveal the origin of 
the disk wind with the best spectroscopic data ever obtained 
for absorption-line diagnostics of the highly ionized plasma in 
GRS 1915+105. Secondly, we study the origin of the contin- 
uum emission in this state by utilizing the ideal combination 
of Chandra and RXTE. The simultaneous HETGS data enable 
us to take into account all the spectral features in the analysis 
of the RXTE PCA spectra, including interstellar/circumstellar 
absorptions with anomalous metal abundances, narrow ab- 
sorption lines by highly ionized gas, and a broad iron-K emis- 
sion line. § 2 summarizes the observations and data reduction. 
The results and discussion are presented in § 3 and § 4, respec- 
tively. We assume the distance toward GRS 1915+105 of 12.5 
kpc and the inclination of 70° throughout the paper. Quoted 
errors for spectral parameters are 90% confidence level for 
a single parameter unless otherwise mentioned. Multiwave- 
lengths properties of GRS 1915+105 including radio and in- 
frared data over a long period covering our observation are 
presented by lArai et all (l2008h . 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra Data 

The Chandra HETGS observation of GRS 1915+105 was 
carried out from 2007 August 14 5:22 (UT) to 19:35 for an 
exposure of 47.4 ksec as a target-of-opportunity (ToO) ob- 
servation (obsID 7485). Figure 1 shows the long term light 
curves of GRS 1 9 1 5+ 1 05 monitored by the RXTE ASM in the 
2-12 keV band, its hardness ratio HR2 between the 5-12 keV 
and 3-5 keV bands, and that by Swift BAT in the 15-50 keV 
band. To catch the source in "soft state", we set the trigger 
condition for the Chandra ToO that the following two crite- 
ria holds for 2 or more successive days: (1) the ASM count 
rate < 30 c/s and (2) HR2 < 1.1. As seen from the figure, 
GRS 1915+105 made a state transition around 2007 August 
8, entering into a quasi-stable "soft state" where the 15-50 
keV flux became below the detection limit of the BAT. The 
position of the arrows denote our observation epoch. To make 
the effects of pile-up least, we adopt graded telemetry mode 
and readout only one-third of the full frame using a subarray. 
We did not use ACIS-S0 and S5, because GRS 1915+105 is 
heavily absorbed and has little photons below — 1 keV. 

We reduce the data using the CIAO version 3.4 and the 
CALDB version 3.4.5 provided by the Chandra X-ray Cen- 
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FIG. 1. — Upper: the long term light curve of GRS 1915+105 observed 
with the RXTE ASM. Middle: the hardness ratio (HR2) of the RXTE ASM 
count rates between 5-12 keV and 3-5 keV bands. Lower: Swift BAT light 
curve in the 15-50 keV band. The epoch of our observations is denoted by 
the arrows. The dashed lines correspond to 30 cts/s (upper) and HR2 =1.1 
(middle), thresholds for triggering our ToO observations aiming at the "soft" 
state. 

ter (CXC). To determine the zero-th order position most 
accurately under heavy pile-up, we utilize the tgjindzo 
script 4 . The accuracy of relative wavelengths calibration of 
the HETGS is 0.05 pixel in position within a chip , corre- 
sponding to 0.00028 in the HEG (or, 14 km s" 1 at 2 keV and 
47 km s" 1 at 7 keV), although the accuracy may be worse if 
we use the data over multiple chips, because of a chip-gap er- 
ror (Marshall, H, private communication). As described be- 
low, we are particularly interested in the energy determination 
in the 2.0-7.0 keV range (between the Ka lines from Si XIV 
and Fe XXVI), which is covered within a single chip in the -1 
order spectrum. We confirm that the +1 data give consistent 
results with the -1 one. Thus, velocity difference between Si 
XIV and Fe XXVI ions can be measured within an accuracy 
of <50kms _I . 

In this paper we mainly utilize the lst-order HEG data, 
which cover energies above 4 keV and has twice better en- 
ergy resolution than the MEG. For spectra below 1.6 keV, 
the MEG data are utilized instead, which have about three 
times larger effective area there. Pile-up is correc ted accord- 
ing to the method described in lUeda et ail d2004l) ; we cal- 
culate the fractional count rate loss as / = 1 - exp(-7.5/?), 
where R is the count rate of the ACIS pixel in units of count 
per pixel per frame, and multiply the observed spectra by 
1 /(l — f) ~ 1+7.5/?. The correction factor is found to be 15% 
at maximum around 4-5 keV in the HEG data. The spectra 
from the +1 and -1 orders are summed. We conservatively 
add a systematic error of 3% in each bin of the spectra. The 
XSPEC package (version 1 1.3.2) is used for spectral analysis. 

2.2. RXTE Data 

Following the trigger of the Chandra ToO, Rossi X-ray 
Timing Explorer (RXTE) performed two ToO observations of 
GRS 1915+105 at 2007 August 14 07:36-08:57 (Epoch F) 
and 17:00-18:03 (Epoch IF) during the Chandra observa- 
tion. RXTE carries two pointing instrumen ts, the Proportional 
Counter Array (PCA: iJahoda et ail 120061) and the High En - 
ergy Timing Experiment (HEXTE: iRothschild et al.l 1998). 

4 available from the MIT/CXC home page, 

http://space.mit.edu/cxc/analysis/findzo/findzo.html 
5 http ://space . mit.edu/CXC/calib/hetgcal .html 
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We only analyze the PCA data in this paper, because the spec- 
tra of GRS 1915+105 were very soft in our observations and 
hence the HEXTE data do not give strong constraints. The 
HEADAS package (version 6.4) is used for the analysis. The 
data are selected by the following criteria: 1) RXTE was not 
in the South Atlantic Anomaly, 2) the elevation angle was be- 
low 10 degrees, and 3) the offset angle between the pointing 
direction and the target was less than 0.02 degree. The net 
exposures are 2.4 ksec and 3.2 ksec for Epochs F and IF, re- 
spectively. 

We utilize the "standard 2" PCA data for spectral analysis, 
which have 16 sec time resolution and 129 energy channels. 
We extract energy spectra from the data of the top layer of 
Proportional Counter Unit (PCU) 2, the best calibrated one. 
The background is subtracted using the model applicable for 
bright sources. The deadtime, about 3^1%, is corrected. A 
systematic error of 1.5% is added for each energy channel. 
To check the calibration accuracy, we also analyze the spectra 
of Crab Nebula observed on 2007 August 11. We confirm 
that the combined PCU 2 and HEXTE (Cluster-B) spectra in 
the 3.5-200 keV band are well reproduced by a single power- 
law with a photon index of 2.07±0.01 for a fixed Galactic 
absorption of yV H = 3 x 10 21 cm -2 . The flux in the 3-20 keV 
band determined by PCU 2 is 2.63 x 10~ 8 ergs cm" 2 s" 1 , which 
is con sistent within 10% of the nominal value dToor & Sewardl 
[1971 . 

For timing analysis of the PCA, we use two types of single- 
bit data covering the energy range of 3.2-5.6 keV and 5.6- 
14.8 keV, and event-mode data covering energies above 14.8 
keV. They have a time resolution of 122 /is (single bit) and 
16 /is (event mode), respectively. The light curves from these 
data are co-added to produce one light curve in the 3.2-37.9 
keV band. The power spectral density (PSD) is then calcu- 
lated by the powspec ftool in the XRONOS package. It is 
known that the PSD computed from PCA data in the Leahy 
normalization is less than 2, a value expected from Poisson 
distribu tion, at high freque ncies above 10 Hz due to deadtime 
effects (Zhang et al. 1995). We thus subtract the dead-time 
corrected Poisson noise level from the observed PSD, f ollow- 
ing the method by Morgan, Remilla rd. & Greinerl (119971) . 

3. RESULTS 

3.1. Light Curves and Power Spectra 

Figure 2 shows the light curves of GRS 1915+105 with 16 
sec resolution obtained with the Chandra HETGS in the 1-8 
keV band (upper panel), RXTE PCA in the 3.2-38 keV band 
(middle), and the hardness ratio between the 5.6-38 keV and 
3.2-5.6 keV PCA count rates (lower). In the upper panel, we 
plot the sum of the MEG and HEG count rates including the 
1st, 2nd, and 3rd orders. As noticed from the figure, the av- 
eraged flux level was roughly constant in the first half and 
then started to increase from the middle of the observation, 
accompanied with spectral hardening. The PCU2 count rate 
increased by about 70% in the latter observation. It is also 
seen that the flux is highly variable on a shorter time scale 
than 10 4 sec. For later discussion, we divide the epoch into 
two, Epoch I (before 2007 August 14 13:06:40) and Epoch 
II (after). We refer the two intervals of the RXTE observa- 
tions as Epoch F (former) and IF (latter). These observations 
would be classified a s "<t>" ( steady State A) in the nomencla- 
ture of lBeilonietail ((2000), since the hardness ratio values 
are low, there is limited variability, and there is no structure 
in the color-color diagram (not shown). During Epoch II, the 
count rate approaches the level that is more typical of the "<5" 
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FIG. 2. — Upper: the Chandra HETGS light curve (sum of the HEG 1st, 
2nd, and 3rd order events) with a time resolution of 16 sec. No correction 
is made for pile-up. Middle: RXTE PCA light curve in the 3.2-38 keV band 
with a time resolution of 16 sec (only PCU2 data are utilized). Lower: the 
PCA hardness ratio between the 5.6-38 keV and 3.2-5.6 keV bands. 
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FIG. 3. — The normalized power spectral density in Epochs I' (upper) 
and IF (lower) obtained from the RXTE PCA data covering the 3.2-38 keV 
band. The Poisson noise level was subtracted after dead time correction. The 
position of the 67 Hz QPO is denoted by the arrow. 

class, where there are transitions between States A and B, and 
hence Epoch II can be considered to be an interval of "bright 
<j>" conditions. 

The power spectral density in the 3-38 keV band calcu- 
lated from the RXTE PCA is shown in Figure 3. While there 
are no clear QPO features in Epoch F, a QPO is detected at 
66.8±1.1 Hz at 4.3er level in Epoch IF. The coherence pa- 
rameter in this QPO is Q = ^o/FWHM « 1 1, where fo is the 
central QPO frequency and the fit uses a Lorentzian profile. 
We note that the difference in the QPO intensity between the 
two epochs is not highly significant, however. The integrated 
fractional root mean square around the 67 Hz QPO is found to 
be 0.8%±0.4% and 1.6%±0.2% in Epoch F and IF, respec- 
tively (the errors are la). 



4 



GRS 1915+105 in soft state 



TABLE 1 

Spectral Parameters of the Time-averaged 
HETGS Spectrum 



Parameter 



Best-fit 



Column Densities" (Nn) 



H(10 22 cm 4 ) 


2.78±0.03 


Mg (10 22 cnT 2 ) 


6.3±1.4 


Si (10 22 cnr 2 ) 


5.8±0.4 


S (10 22 cm" 2 ) 


7.6±1.2 


Fe (10 22 cnT 2 ) b 


10±1 


Continuum 


A 


7.15±0.25 


r 


1.17±0.06 


£foid (keV) 


3.80±0.17 


Iron-K Emission Line d 


E cm (keV) 


6.55±0.04 


lcr width (keV) 


0.17±0.04 


E.W. (eV) 


45±8 


£cen (keV) 


6.54±0.03 



lcr width (keV) 
r m (r g ) 
rout (>g) 

P 



0.05 (fixed) 

600+ 300 
ouu -200 

100000 (fixed) 
-3 (fixed) 



NOTE. — The errors are 90% confidence level for a single 
parameter. 

a Equivalent hydro gen column d ensities by assuming the solar 
abundances I Anders & Grevesse 1989) between the element 
and hydrogen (the abundance ratios within each group of H- 
He-C-N-O, Ne-Na-Mg-Al, S-Cl-Ar-Ca, and Cr-Fe-Co-Ni are 
fixed at the solar values). The columns of Mg, Si and S are 
determined by the spectral features around each K edge, while 
that of H is estimated from the fit to the whole HEG first or- 
der spectrum by taking account of dust scattering based on 
the Draine (2003) cross section (see text). b Determined from 
the Suzaku data. A cutoff power law model with the form of 
A E~ T exp(-£ / '£f id) is adopted, where A is the normalization 
at 1 keV in units of photons cm -2 s keV -1 . d Modeled either 
by a single gaussian (upper section) or by a gaussian blurred 
with the diskline profile (lower). 



3.2. HETGS Spectra 

Figure 4 shows the HETGS spectra integrated over the 
whole observation. Here we only correct the spectrum for 
effective area as a function of energy, but not for the instru- 
mental energy resolution, which is 3.2 eV at 2.0 keV and 40 
eV at 7.0 keV in the HEG. The MEG data are utilized in the 
1.0-1.6 keV band, while the HEG data for the rest. Although 
the overall spectral shape changed from Epoch I to II as men- 
tioned above, we first study this spectrum as the base data 
from which we can best constrain local spectral features in 
time-averaged sense. We find that the overall continuum can 
be well described by a cutoff power law model, expressed as 
A E~ T exp(-E I £foid), subject to heavy interstellar absorption 
with non-solar abundances. The best-fit parameters are sum- 
marized in Table 1 . Note that the cutoff power law model is 
only empirical and does not have physical meaning itself; we 
will model the continuum emission by physical models for the 
PCA spectra in §3.3. 

It is a critical issue to accurately determine the column den- 
sities of major elements in the interstellar (plus circumstellar) 
matter toward GRS 1915+105 for detailed study of its spec- 
tral model. In spectra with insufficient energy resolution, a 
deep absorption edge produced by interstellar matter can be 



coupled with other spectral features such as emission lines at 
similar energies. It could also affect the continuum shape at 
low energies where the photo-electronic absorption cross sec- 
tion is large. The Chandra high resolution data provide us 
with the best opportunity to solve this issue for this source. 
The result also has an important implication for study of ele- 
mental abundances in th e Galaxy, although we do not discuss 
it here. iLee et al] d2002l) determined column densities of Mg, 
Si, S, and Fe by measuring the K-edge depth of each element, 
which are quite useful a nd are adopted b y some authors for 
spectral modelling (e.g., Done et al. 20041). 
In this paper, we re-determine those val ues from o ur dat a 



that have a better statistics than that o f ILee et al] ( 2002 ) 
following the same method described in lUeda et all (|2005). 



We adopt the ph oto-electric absorption cross section by 
Wilm s~etaT] (l2000t) . which is available as TBvarabs model in 
XSPEC. The column density of lighter elements than Mg is 
determined by a spectral fit to the whole continuum. Here 
we consid er the effect of dust scattering based on the cross 
section by Draine (2003). Those of other elements are lo- 
cally fit around the K-edge by taking into account X-ray ab- 
sorption fine structure (XAFS); in particular, it is difficult to 
accurately estimate the column density of Si just by measur- 
ing a depth of Si-K edge, due to the complicated XAFS that 
is most likely attributed to silicates grains. We thus model 
it using an exper imental cross section of Si02 as done in 
lUeda et al.1 (I2005I) 6 . Table 1 lists the Hydrogen-equivalent 
column d ensities of each element assuming the solar abun- 
dances by lAnders & Grevessd(ll989l) for each element, which 
are used in our subsequent analysis at their best-fit values. 
The relative abundances within the groups of H-He-C-N-O, 
Ne-Na-Mg-Al, S-Cl-Ar-Ca, and Cr-Fe-Co-Ni are fixed at the 
solar values. We fin d the column de nsity of Si is not as large 
as that estimated by lLee et alJ (120021) . For the column density 
of Fe, we adopt a value determined with Suzaku (Ueda, Y. et 
al. in preparation) because the Chandra result has a large sta- 
tistical error, yet consistent with the Suzaku result, due to the 
limited energy coverage and area above 7.1 keV. 

As seen from Figure 4, we detect 33 narrow absorption lines 
over the whole energy range of HETGS, for which we sim- 
ply model with negative gaussians in the plot. We find that 
32 of them correspond to transition in highly ionized (H-like 
or He-like) ions. The absorption line at 2.470 keV, identi- 
fied as S II, is most probably produced by interstellar matter 
in the Galaxy, which is commonl y seen in the h igh resolu- 
tion spectra of Galactic sources dUeda et alj|2005l) . The en- 
ergy, lcr width, and equivalent widths are listed in Table 2 
with its identification. For weak lines whose width cannot 
be well determined, we assume a velocity dispersion of 100 
km s" 1 or 200 km s -1 , depending on element. In addition to 
those narrow absorption lines, a broad emission line feature 
centered at 6.55±0.04 keV is found, attributable to Ka lines 
of mildly ionized Fe ions. This can be modeled by a sin- 
gle broad (lcr width of 0. 17±0.04 keV) gaussian or by a nar- 
rower (0.05 keV) ga ussian blurred with the "disk line" profile 
(Fabia n et alj|1989l) . which is a more physical interpretation 
for an iron-K line from the accretion disk (see § 4.2 for dis- 
cussion). In the disk-line fit, we assume an outer radius of 
10 5 r g with an emissivity law of r@ where r is the radius and 
(3 = -3, and make the center energy and inner radius as free 
parameters. The parameters of the emission line are summa- 

6 empirically adding two negative gaussians at 1.846 keV and 1.865 keV 
to a single absorption edge at 1.839 keV works as a good approximation. 
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FIG. 4. — The Chandra HETGS spectrum integrated over the whole observation (corrected for the effective area but not for the energy resolution). From upper 
to lower: 1.0-1.6 keV, 1.6-2.2 keV, 2.2-2.8 keV, 2.8-3.4 keV, 3.4-5.6 keV, 5.6-7.8 keV. The HEG lst-order events are used except below 1.6 keV, where the 
MEG 1 st-order data are utilized. The solid curve (red) represents the best-fit continuum (modelled by an absorbed cutoff power law) modified by absorption lines 
(negative gaussians) and a broad iron-K emission (a positive gaussian). The upper axis of each panel marks wavelengths in units of A (blue). 



rized in Table 1 . 

To determine the physical parameters of the ionized plasma, 
we next fit the absorption line profile of major elements (Si, 
S, Ar, Ca, Cr, Mn, Fe) by a Voigt function in the local en- 
ergy band. The Voigt profile, whose opacity is expressed by 
a convolution of a gaussian and a lorenzian (natural broad- 
ening), is a physically more correct model rather than a neg- 
ative gaussian. To do this, we utilize a local m odel (called 
"Kabs ") implemented on XSPEC as described in lUeda et ail 
(2004), where its full formula is given. The free parame- 
ters are three per each ion, a column density (cm -2 ), a line- 
of-sight velocity dispersion b (km s" 1 ), and a Doppler shift 



z = v/c, where v is the mean bulk velocity of plasma (nega- 
tive for outflow) and c is the light velocity. The atomic data 
are taken from the database provided by National Institute of 
Standards and Technology (NIST) 7 . In this paper, we treat 
an ideal case where contribution of re-emission line from the 
plasma ou t of the line - of-s ight is neglected, corresponding to 
Case I in Ueda et al. (2004) or "non-radiative de-excitation 
limit" in iMasai & Ish ida (2004J). It is because the contribu- 
tion of the emission cannot be well constrained by our data, 
being consistent with zero within the statistical errors; when 

7 Version 3.1.5, http://physics.nist.gov/asd3 [2008 September 6] 
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GRS 1915+105 in soft state 



TABLE 2 

List of Absorption Lines Observed in the HETGS spectra 



£obs 


1<t width 3 


E.W. 


line ID b 


£theoiy (term of final state) c 


(keV) 


(eV) 


(eV) 




(keV) 


1.277(4) 


0.43* 


1.8 


NeX4p 


1.27707 ( 2 P 1/2 ), 1.27713 ( 2 P 3 / 2 ) 


1.4734(4) 0.001« 1.0) 


1.7 


Mg XII 2p 


1.47169 ( 2 P 1/2 ), 1.47264 ( 2 P 3/2 ) 


L.DJZ.Z{0) 


l.U(OJ 


1.4 


FeXXIV 2s^4p 1.55091 ( 2 P 1/2 ), 1.55289 ( 2 P 3 / 2 ) 


I Tin i p~i\ 


Hf 1 \ 

2(1) 


1.9 


Al XIII 2p 


1.72769 ( 2 P, n ), 1.72899 ( 2 P 3/2 ) 








Fe XXIV 2s^5p 1.72876 ( 2 Pi/ 2 ), 1.72978 ( 2 P 3/2 ) 




A A/I / ^ 1 1 \ 


0.67 


Mg XII 3p 


1.74456 ( 2 P 1/2 ), 1.74484 ( 2 P 3 / 2 ) 


l.oLo{Z) 


U.Ol 


0.31 Ni XXVI 2s-*4p 1.81758 ( 2 P 1/2 ), 1.82033 ( 2 P 3 / 2 ) 


l.oZ /op J 


U.Ol 


0.75 Fe XXIV 2s^6p 1.82640 ( 2 P l/2 ), 1.82698 ( 2 P 3 / 2 ) 


1 .042(2) 


U.01^ 


0.54 


Mg XII 4p 


1.84002 ( 2 P 1/2 ), 1.84014 ( 2 P 3/2 ) 


l.oool(o ) 


0.63* 


0.49 


Mg XII 5p 


1.88419 ( 2 P,/ 2 ), 1.88426 ( 2 P 3/2 ) 


2.0063(3) 


1.4(2) 


3.4 


Si XIV 2p 


2.00433 ( 2 Pi /2 ), 2.00608 ( 2 P 3/2 ) 


2.029(2) 


0.68* 


0.39 Ni XXVI 2s-+5p 2.02775 ( 2 Pi/ 2 ), 2.02916 ( 2 P 3/2 ) 


2.3777(5) 


0.7(< 1.8) 


1.7 


Si XIV 3p 


2.37611 ( 2 P r/2 ), 2.37663 ( 2 P 3/2 ) 


2.460(4) 


0.82* 


0.54 


SXV 2p 


2.46063 ('Pi) 


O /I *7 ft / 1 \ 

2.4/0(1) 


2(1) 


1.3 


S II 3p d 


2.4694 d 


Z.JU /( 1 J 


A ST* 


1.2 


SiXTV4p 


2.50616 ( 2 P 1/2 ), 2.50638 ( 2 P 3/2 ) 


2.568(3) 


0.85* 


0.46 


Si XIV 5p 


2.56632 ( 2 Pi /2 ), 2.56644 ( 2 P 3/2 ) 


2.599(4) 


0.87* 


0.44 


Si XIV 6p 


2.59899 ( 2 Pi/ 2 ), 2.59906 ( 2 P 3/2 ) 


2.6229(3) 


2.5(3) 


5.2 


S XVI 2p 


2.61970 ( 2 Pi/ 2 ), 2.62270 ( 2 P 3/2 ) 


3.1077(7) 




2.0 


S XVI 3p 


3.10586 ( 2 Pi /2 ), 3.10675 ( 2 P 3/2 ) 


3.277(2) 


1.1* 


0.92 


S XVI 4p 


3.27589 ( 2 Pi/ 2 ), 3.27627 ( 2 P 3/2 ) 


3.3227(6) 


2.9(7) 


4.5 


Ar XVIII 2p 


3.31818 ( 2 Pi/ 2 ), 3.32299 ( 2 P 3/2 ) 


3.355(4) 


1.1* 


0.44 


S XVI 5p 


3.35454 ( 2 P 1/2 ), 3.35473 ( 2 P 3/2 ) 


3.905(4) 


1.3* 


1.1 


Ca XIX 2p 


3.90226 ('Pi) 


3.937(4) 


1.3* 


1.2 


Ar XVIII 3p 


3.93429 ( 2 P 1/2 ), 3.93572 ( 2 P 3/2 ) 


4.1070(8) 


5(2) 


6.1 


Ca XX 2p 


4.10015 ( 2 P 1/2 ), 4.10750 ( 2 P 3/2 ) 


4.151(6) 


1.4* 


0.80 


Ar XVIII 4p 


4.14974 ( 2 P 1/2 ), 4.15034 ( 2 P 3/2 ) 


4.865(5) 


1.6* 


1.7 


Ca XX 3p 


4.86192 ( 2 P 1/2 ), 4.86410 ( 2 P 3/2 ) 


5.684(6) 




3.3 


Cr XXIII 2p 


5.68205 ('Pi) 


5.931(4) 


7« 16) 


5.8 


Cr XXIV 2p 


5.91650 ( 2 P 1/2 ), 5.93185 ( 2 P 3/2 ) 


6.17(2) 


4 i** 


3.6 


Mn XXIV 2p 


6.18044 ('P[) 


6.45(1) 


nO 


7.9 


Mn XXV 2p 


6.42356 ( 2 P 1/2 ), 6.44166 ( 2 P 3/2 ) 


6.692(3) 


28(3) 


40 


Fe XXV 2p 
Fe XXIV 2p 


6.70041 ('P) 
6.67644 ( 2 P 1/2 ), 6.67915 ( 2 P 3/2 ) 


6.975(2) 


16(2) 


40 


Fe XXVI 2p 
Cr XXIV 3p 


6.95196 ( 2 P 1/2 ), 6.97317 ( 2 P 3/2 ) 
7.01726 ( 2 P 1/2 ), 7.02181 ( 2 P 3/2 ) 



NOTE. — The single number in the parenthesis indicates the statistical error in 
the last digit (90% confidence level for a single parameter), 
a ^ **. fj xe( j at a va i ue coiresponding to a velocity dispersion of 100 km/s (*) 
and 200 km/s (**). b The initial state is "Is" unless otherwise indicated. Taken 
from the NIST database (version 3.1.5) except for the S II 3p line, for which 
an observed value by Ueda et al. (2005) is shown. d Most probably due to the 
interstellar gas (see text). 

we include such emission lines emitted at r = (1 — 3) X 10 5 r g 
with an emissivity law of r" 3 (corresponding to the case of an 
optically thin, self-similar flow with a constant velocity) and 
an inclination of 70° in the spectral model, we obtain an upper 
limit for the solid angle of the plasma as H/2tt < 0.1 1 (90% 
confidence level) from a simultaneous fit to the Si XIV and S 
XVI absorption lines. We find that inclusion of these emis- 
sion lines at the maximum level could increase the column 
density obtained for the same velocity dispersion by «30% 
for Si XIV and ss50% for Fe XXVI, which are within the sta- 
tistical errors. 

We hereafter analyze the HETGS spectra of Epochs I and II 
separately, since the absorption line profile could change ac- 
cording to the evolution of the continuum spectrum. Figure 5 
shows the HEG lst-order spectra corrected for effective area 
(i.e., incident spectra except for blurring with the energy res- 
olution). It is evident that the spectrum is harder in Epoch II 




1 2 5 



Energy (keV) 




Energy (keV) 



FIG. 5. — Upper: The HEG lst-order spectrum (corrected for effective 
area) in Epoch I. Lower: The HEG 1 st-order spectrum (corrected for effective 
area) in Epoch II. 

than in Epoch I. In this subsection, we focus on the absorption 
lines. Detailed analysis of the continuum is given in next sub- 
section using the RXTE PCA data, which have a wider energy 
coverage than the HETGS. 

To best constrain the physical parameters of each ion, we 
fit two absorption lines of Ka and K(3, whenever available 
(i.e., except for Mn and Fe), simultaneously in the local 
bands. This works quite well to de-couple the inevitable anti- 
correlation between a column density and a velocity disper- 
sion in the curve-of-growth analysis. We simply adopt an ab- 
sorbed cut-off power model for the continuum by fixing 
and the interstellar absorption. The velocity dispersion in two 
epochs is found to be consistent with a constant value within 
the errors. We thus assume the best-fit velocity dispersion ob- 
tained from the summed spectrum of Epochs I and II when 
deriving the column density and Doppler shift of the Si, S, 
Ar, and Ca ions in each epoch. In the iron-K band, we simul- 
taneously fit Ka and K(3 lines of Cr XXIV, Ka of Mn XXV, 
and Ka of Fe XXVI by assuming a same velocity dispersion 
and Doppler shift. The Cr XXIV Kf3 line (7.02 keV) contam- 
inates a high energy part of the Fe XXVI Ka feature in the 
HETGS spectra. We treat the Fe XXV absorption line inde- 
pendently, considering possible uncertainties in the incident 
line energies 8 and a contamination of less ionized ions such 
as Fe XXIV, although we assume the same velocity disper- 

8 Unlike a single electron system in H-like ions, there could be more un- 
certainties in the atomic database for He-like (and more electrons) ions. 
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sion as that of Cr XXIV, Mn XXV, and Fe XXVI. We note 
that the results of Fe XXV may not be reliable to be taken 
at their face values but just provide a rough estimate. For Si 
XIV and Fe XXVI, we introduce an upper limit for the col- 
umn density constrained by the corresponding K-edge depth 
measured by HETGS and by PCA, respectively; otherwise, 
the fit to the absorption line profile leads to unphysically large 
column densities. 

Figure 6 shows a blow-up of absorption line profile in 
Epochs I and II for each set of elements around their Ka lines. 
The solid curves represent the best-fit model. Table 3 gives the 
summary of the parameters. We find that all the ions have sig- 
nificant blue shifts, w 150 km s" 1 for Si XIV and «500 km s" 1 
for Fe XXVI. The difference of the velocity between H-like Si 
(and S, Ar, Ca) and Fe ions is robust, larger than the system- 
atic error in the relative wavelength calibration as mentioned 
in § 2.1. The velocity dispersion increases with the atomic 
number, from 70 km s" 1 for Si to w 200 km s" 1 for Fe. From 
Epoch I to II, the column density of all the ions decreases, 
indicating that the plasma was in a higher ionization stage in 
the latter epoch than the former; the number ratio between Fe 
XXV and Fe XXVI is also consistent with this picture. 

3.3. Continuum Model 

Once the local spectral features and interstellar absorption 
are accurately determined with the Chandra data, we can now 
study the continuum emission with the RXTE PCA spectra in 
the 3-25 keV range with least uncertainties even if they can- 
not be resolved with the PCA data alone. We assume that the 
PCA spectra integrated in Epochs I' and IF well represent the 
HETGS spectra in Epochs I and II, respectively, although the 
time coverage by RXTE is shorter than that by Chandra. In 
the analysis of PCA spectra, we fix the parameters of absorp- 
tion/emission lines and absorption column densities at the the 
best-fit values obtained above. We introduce four additional 
spectral features that reside outside the HETGS band: absorp- 
tion lines at 7.84 keV and 8.20 keV and edge structures at 
8.6 keV and 9.0 keV. The absorption lines correspond to Fe 
XXV K/3 + Ni XXVII Ka and Fe XXVI K/3 + Ni XXVIII 
Ka, for which we a ssume an equivale nt width of 30 eV and 
20 eV, respectively dKotani et al.ll2000l) . The edge models ap- 
proximately represent spectral modification by the bound-free 
(photo-absorption) and bound-bound transition close to the K- 
edge energy of Fe XXV and Fe XXVI ions, respectively. The 
ratio of the edge depth between Fe XXV and Fe XXVI is fixed 
at 1.1 (Epoch F) and 0.44 (Epoch IF), based on the results 
from the HETGS. In the continuum modeling, we always in- 
clude a reflection component that is expected from the broad 
iron-K emission line det ected with the HETGS. We model 
it with the pexriv code dMagdziarz & Zd ziarski 1995), or its 
convolution model. We fix the inclination, temperature, abun- 
dances of Fe, and that of other metals to be 70°, 10 6 K, 3.6 
solar, and 2.1 solar, respectively. To make it consistent with 
the central energy and equivalent width of the iron-K emis- 
sion line observed with the HETGS, the ionization parameter 
£ re fi is fixed at 20 (Epoch F) and 40 (Epoch IF), and the solid 
angle of the reflector f2 re fl/27r is restricted be less than 0.5. 
Since we consider the accretion disk as the reflector, it is fur- 
ther blurred with the diskline kernel assuming an emissivity 
profile of r" 3 from radii of (600- 10 5 )r g , again based on the 
results of the HETGS for the iron-K line. 

We find that a cutoff power law model give a good (math- 
ematical) description for the overall PCA spectra in the 3-25 
keV range as is for the HETGS spectra. As a physical model, 



we fir st try the multi color disk (MCD) model ( Mitsu da et al.l 
1984) plus a power law, which has been widely used to fit 
the spectr a of canonical black holes in the high/soft state 
(see e.g., iTanaka & Lewin|[l995l) . This model is not accept- 
able for our PCA spectra of GRS 1915+105, with reduced 
X 2 (d.o.f) of 2.19 (44) and 3.98 (44) for Epoch F and IF, re- 
spectively. Note that this model can however fit the HETGS 
spectra whose coverage is limited below 7.5 keV, demonstrat- 
ing that continuum modeling only by using data below <10 
keV are sometimes unreliable. Next, we try a power law 
plus a "p-free disk" model, where the temperature profile is 
express as a T = T m {r/r m )~ p with p being a free parameter 
instead of a fixed value of 0.75 in the MCD model. This 
model has been successfully applied for bl ack holes accret- 
ing a t high fractional Eddington r atio (e.g., Mineshig e et al.l 
ll994tlKubota & Makishimall2004l) . We find that the fit is sig- 
nificantly improved (and is even successful for the Epoch IF 
data); we obtain p = 0.43^;|, T m = 2.31 ±0.12 keV, r m = 
5.7^;^ km with reduced x 2 (d.o.f) of 1.46 (43) for Epoch F, 
and p = 0.49^;|, T in = 2.48 ± 0.09 keV, and r in = 6.1^ km 
with reduced \ 2 (d.o.f) of 1.01 (43) for Epoch IF. Here r; n is 
corrected for cosine of the inclination by assuming the emis- 
sion region is more like a disk rather than spherical geometry. 

We finally try a Comptonization model, which has been 
used to mo del the spectra of black hol es in the very high 
state (e.g.. iKubota & Makishima 2004) and the low/hard 
state (e.g., iGirlinski et al.1 1 19971) . This model is found to 
be a good des cription of the sp ectra of GRS 1915+105 in 
severa l cases (|Done et al.l l2004t). Here we utilize thComp 
code (IZycki. Do ne. & Smith 1999), which is appropriate even 
when the optical depth of Comptonizing corona is larger than 
3, as in our case. The spectrum of seed photons is assumed to 
be the MCD model. The free parameters are the temperature 
of Comptonizing electrons T e , the photon index T, with the in- 
nermost temperature T m and radius r m for the incident MCD 
component. We can infer r; n from the conservation of photon 
numbers in the Comptonization process (see § 4.2). The scat- 
tering optical depth can b e calculated from T R and T ac cording 
to the following formula ( iSunvaev & TitarchukHl980l) : 



r e y 2 - 25 + (7;/511keV) . [(r + 05) 2_2.25] 5 ' 

Since we find no direct (i.e., without Comptonization) MCD 
component, we do not included it in our mo del, unlike the 
ca se of GRO J 1655-40 dTakahashi et aIT2 008) and Cygnus X- 
1 (Ma kishima et al.l l2008). The upper limit on the fraction of 
the non-Comptonized flux to the Comptonized one (the flux 
before being Comptonized) is 20% (Epoch F) and 3% (Epoch 
II). 

The fit with the Comptonization model is found to be much 
better than the p-free disk model, yielding reduced \ 2 (d.o.f) 
of 0.94 (43) and 0.56 (43) for Epoch F and IF, respectively. 
Figure 7 (upper) show the PCA spectra folded with the energy 
response, with the best-fit model (solid line) and residuals in 
units of x- Figure 7 (lower) show the incident photon spectra 
in the form of E 2 Fe, where E and Fe is the energy and pho- 
ton spectrum, respectively. The contribution of the reflection 
component including an iron-K emission line calculated self- 
consistently is separately plotted. The best-fit parameters (in- 
cluding those of the reflection component) are summarized in 
Table 4. We find the Comptonizing corona has a low electron 
temperature of «3.6 keV and a large scattering optical-depth 
of m 5, which are constant in the two epochs within the error. 
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Epoch II 




3.2 
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channel energy (keV) 
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Epoch I 



4.05 4.1 4.15 4.2 4.25 4.3 3.9 3.95 

channel energy (keV) 
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Photons/cm : 
0.1 
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FIG. 6. — A blowup of the time-averaged HEG lst-order spectra around each K-absorption line in Epochs I (left) and II (right). The best-fit model employing 
Voigt functions is plotted with a line (red). From upper to lower: Si XIV, S XVI, Ar XVIII, Ca XX, Cr XXIV + Mn XXV + Fe XXV + Fe XXVI. 



TABLE 3 

Parameters of the Ka Absorption Lines Fit by Voigt Profile 



Epoch I 



Epoch II 



Ions 



(xi<r 3 ) 



(kms" 1 ) 



(10 18 



cm ) 



(xlfr 3 ) 



b b 
(km s~ 



AV, 



) 



(10 11 



Fe XXVI 
Fe XXV C 
MnXXV 
Cr XXIV 
CaXX 
Ar XVIII 
SXVI 
Si XIV 



1.68 ±0.18 
0.66 ±0.07 
(= Fe XXVI) 
(= Fe XXVI) 
061 +o.i3 
u - -0.19 
0.41+0.10 

51+^ 
-Q.Q3 

u - JJ -0.04 



190+|> 
(= Fe XXVI) 
(= Fe XXVI) 
(= Fe XXVI) 
110+30°* 

105+** 
74+'°* 



13+ 



3+7 

-8 

7.2±1.4 
0.19±0.06 
0.19±0.04 

31+012 

0.26±0.05 

47+0.08 
"■"-If 
0.47+06 



1.72+ 



7+0.09 
-0.05 

0.31 ±0.08 

(= Fe XXVI) 

(= Fe XXVI) 

63+°- 17 

62+ 20 
-Q15 
55+0- 08 

0.50 ±0.06 



210±10 
(= Fe XXVI) 
(= Fe XXVI) 
(= Fe XXVI) 
(= Epoch I)* 
(= Epoch I)* 
(= Epoch I)* 
(= Epoch I)* 



10+0 
2.2±0.3 
0.16±0.06 
0.12±0.05 

n 1 0+0.04 

0.14±0.02 

23+ - 07 
-QQ3 
35+0 ™ 
u ■ J -0.08 



NOTE. — The errors are 90% confidence level for a single parameter. 
a Blue shifts (i.e., outflow velocity relative to the light velocity). b Line-of-sight velocity dispersion. There 
may be systematic uncertainties in the parameters of Fe XXV (see text). 
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TABLE 4 

Best-Fit Parameters of the RXTE PCA Spectra 



Parameter 


Unit 


Epoch P 


Epoch IF 


'"in 
Tin 
T 

r 

n refl /27r b 


(keV) 
(Kev ; 

(erg cm s" 1 ) 


1 13±7 
0.87±0.06 

J - 1S -0.15 

4.8±0.3 

5^>-° 

u - J -0.06 
20 (fixed) 


79±8 
1.11 ±0.06 

2.99±0.14 
5.2±0.4 ) 

42+0.O8 

40 (fixed) 


r cdgc (9.0 keV) c 
t c(1) , c (8.6 keV) c 




0.25±0.02 
0.28±0.03 


0.16±0.03 
0.07±0.02 


F 2 -io d 
£o.oi-ioo e 


(erg cm -2 s -1 ) 
(ergs -1 ) 


1.0 X 10- s 
6.6 x 10 38 


1.6 x 10" 8 
8.3 x 10 38 


X 2 (d.o.f)' 




40.4 (43) 


24.1 (43) 



NOTE. — The errors are 90% confidence level for a single pa- 
rameter. 

a The optical depth of the Compton cloud calculated from 7i and F 
(see text). b A reflection component is included with a solid angle 
f2 lcn and an ionization parameter £ re fl, blurred by the diskline ker- 
nel assuming the emissivity profile of r~ 3 between r m = 600 r g and 
'but = 10 5 r g . c Optical depth of the edge model at 9.0 keV and 8.6 
keV. The ratio between 9.0 and 8.6 keV is fixed at 1.1 (Epoch F) 
and 0.44 (Epoch IF). d Observed flux in the 2-10 keV band. 6 In- 
trinsic luminosity of the direct continuum (without the reflection 
component) in the 0.1-100 keV band, corrected for absorption. 
The Comptonized radiation is assumed to be isotropic/ The \ 2 
value and degree of freedom. 

The hardening in Epoch II can be attributed to the intrinsic 
change of the seed-photon spectrum. The 0.01-100 keV in- 
trinsic luminosity corrected for absorption and reflection is 



found to be 6.6 x 10 38 erg s" 1 (Epoch I') and 8.3 x 10 38 erg 
s" 1 (Epoch II'), by assuming isotropic emission. These cor- 
respond to 0.3-0.4 times the Eddington limit of 14 M black 
hole, and could become 0.9-1.2 times the Eddington if the 
emission region is optically thick and has a disk-like geome- 
try. We discuss the structure of the accretion disk in § 4.2. 

4. DISCUSSION 
4.1. Nature of Accretion Disk Wind 

We have detected at least 32 narrow absorption lines from 
a highly ionized plasma (H-like or He-like) in the HETGS 
spectra of GRS 1915+105 in "soft state" that are the deep- 
est among those detected from the source, achieving our pri- 
mary aim with this ToO observation. This enables us to in- 
vestigate the physical properties of the plasma with the best 
accuracy ever achieved. The velocity and line-of-sight veloc- 
ity dispersion are determined to be 150-500 km s" 1 and 70- 
200 km s , respectively, which differ by element. Since the 
proper motion of GRS 1915+105 is estimated at 3 ± 10 km 
s" 1 from lGreiner. Cuby. & McGaughreanl(l2001h (considering 
the fact that the black hole mass is about 12 times larger than 
the companion in this system), it is revealed that the plasma 
is definitely in outflow, and hence the origin is a wind (most 
likely an accreti on disk wind) as suggested in previous studies 
dLee et alj2 002). The depth of these absorption lines changed 
according to the X-ray flux, consistent with the picture that it 
is photo-ionized by the radiation. 

Estimating the location of the wind is crucial to discrimi- 
nate its origin. We estimate it b y a simple argume nt using the 
ionization parameter £ = L/nr 2 dTarter et"aT1 ll969). where L is 
the luminosity, n is the density, and r is the distance from the 
source. As discussed below, £ may not always constant along 
the line of sight but here we discuss its averaged value for first 



order estimate. For this purpose, we perform simulation of a 
photoionized plasma by the XSTAR code for various values 
of £, assuming our best-fit continuum models in Table 4 and 
the metal abundances in Table 1 . The density n is set to be 
10 12 cm" 3 . We find that the observed number ratio between 
Fe XXVI and Fe XXV ions, 1.8±1.0 and 4.5±1.5, indicates 
log £ = 4.3^6 and log £ = 4.2 ±0.1 for Epoch I and II, re- 
spectively; here L is defined in the 0.0136-13.6 keV band in 
each spectrum (this is the reason why the £ value is larger in 
Epoch I in spite of the lower ionization state). Assuming that 
the length of the wind, Ar, is a similar order to the distance 
from the source, we find nr ~ nAr ~ 10 23 cm" 2 after correct- 
ing for the fraction of Fe XXVII ions and Fe abundance of 3.6 
times the solar. This yields r ~ L/£/nr ~ (2-6) x 10 11 cm, or 
- (1 - 3) x 10 5 r g for a black hole mass of 14 M 9 . This 
well exceeds a possible minimum launching point of ther- 
mally driven wind from an X-ray irradiated disk theoretically 
predicted. The critical radius where the thermal energy of the 
X-ray heated disk exceeds the gravitational energy is given by 



10 I0 (- 



10 s K M 



Tc ^ 

Here Tc is the Compton temperature, at which the net energy 
transfer between electrons and photons is balanced out. In 
our case Tc is estimated to be ~0.3 keV by taking the aver- 
age of photon energy of the continuum in the 0.01-100 keV 
band, and hence r c ~ 4 x 10 12 cm. It is predicted that such 
wind can be launched fr om ~ 0.1 r c or even much smaller 
radii dWoods et alJ fl996). As the luminosity (> 6.6 x 10 38 
erg s" 1 in the 0.1-100 keV band) is about 0.3 times the Ed- 
dington limit, the effect of radiation should be also important 
even if the gas is highly ionized and the opacity by UV-lines 
is little. Thus, the wind we observe is consistent with a ther- 
mally (plus radiation) driven wind and there is no necessity to 
invoke a magnetically driven wind. 

The mass outflow rate carried by the accretion disk wind is 
evaluated as 



M oul = fi„ind'" 2 n'Mpl'„ind 

-(0.6-1.3) x 10 19 ( 



^wind/47T 

^wind 



0.2 



)( 



500 km s- 



-) gs 



Here we have used r 2 n = L/£~(3-6)x 10 34 , m p is the pro- 
ton mass, i7 w ind is the solid angle of the wind, and v w j n d is the 
outflow velocity (which is larger than the observed outflow 
velocity by a factor of cosine of the angle between the direc- 
tion of the wind and our line-of-sight). Thus, the mass outflow 
rate carried by the wind could be comparable to the mass ac- 
cretion rate in the inner part of the disk, which is estimated to 
be (0.7-0.9) x 10 19 g s" 1 (spherical case) or (2-3) x 10 19 g 
s" 1 (disk-like case) by assuming the mass-to-energy conver- 
sion efficiency of 0.1. This demonstrates the importance of 
the wind for understanding the dynamics of the whole disk 
system. A similar situation i s found in the n eutron star low 
mass X-ray binary GX 13+1 dUeda et alj|2004 except for the 
difference of the mass scale. 

We find that H-like ions of heavier elements have higher 
outflow velocity and larger velocity dispersion. In addition, 
there is an implication that Fe XXV shows a smaller velocity 
than that of Fe XXVI, although we have to keep in mind the 

9 Note that here we have assumed that the gas is smoothly distributed i.e., 
that it is not clumped. Clumping would reduce the radius, and hence this 
estimate is an upper limit. 
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o Epoch I 



o Epoch IV 




FIG. 7. — Upper: The RXTE PCA spectra (with error bars) folded with the energy response in Epoch F (left) and IF (right). The upper curve (black) represents 
the best-fit model including the contribution of the reflection component (lower curve, blue). The lower panel shows the residuals of the data from the best-fit 
model, normalized by the error (i.e., in units of x). Lower: the best-fit model in the form of E 2 Fg, where E and Fe is the energy and photon spectrum, respectively, 
in Epoch F (left) and Epoch II" (right). The lower curve (blue) represents the reflection component. 



possible systematic uncertainties in the obtained parameters 
for Fe XXV. These results indicate that we trace different parts 
of the wind by different ions; there is a component of the flow 
that is more highly ionized (i.e., traced by Fe XXVI) and that 
is moving approximately 3 times faster than the gas traced by 
other lines. A ccording to a thermally driven wind model (see 
iNetzer Il2006l and references therein), the velocity of outflow 
can be represented in a power law form of distance, 

v oc r 7 , 

where 7 is in the range of 0-0.5. Then, from the conservation 

of the mass outflow rate, the density is given as 

-2-1 
n oc r ' , 

and therefore 

£ oc r 7 . 

For 7 > 0, ionization becomes stronger at outer parts of the 
wind. In this interpretation, H-like Si ions we observe must 
be located closer to the center where £ is smaller, and hence 
have a smaller outflow velocity, compared with those of H- 
like Fe ions. The larger velocity dispersion found in heavier 
element of H-like ions can be attributed to a larger line-of- 
sight velocity gradient when integrated in a more outer region 
of the wind, probably accompanied by larger turbulent motion 
as well. Constructing a fully quantitative model for the wind 
of GRS 1915+105 is beyond the scope of this paper and is left 
for future studies. 



4.2. Origin of X-ray Continuum Emission in "Soft State " 

The ideal combination of the Chandra HETGS and RXTE 
PCA has enabled us to constrain the continuum model with 
least uncertainties, by correctly taking into account the inter- 
stellar absorption with the (updated) anomalous abundances 
and local spectral features including absorption- and emission 
lines. The spectra of GRS 1915+105 in "soft s tate" (steady 
State A in the definition of iBelloni et al.ll2000h in the 1-25 
keV band can be well explained with a strongly Comptonized 
continuum and its reflection component by an outer part of 
the disk, which is absorbed by a highly ionized disk wind dis- 
cussed above. 

A relativisticall y broadened iron-K emi ssion line as previ- 
ously reported by Marto cchia et al.l d2002l) is not required in 
our data, although it is difficult to perfectly exclude its pres- 
ence from the spectral fit alone, which strongly depends on 
the continuum model (and also on the interstellar absorption). 
The large inner radius we obtained from the line profile with 
the Chandra HETGS, 400r g , does not necessarily mean that 
the accretion disk is truncated at that radii. Rather, as dis- 
cussed below, due to the strong Comptonization by corona 
that completely covers the inner disk, any iron-K features aris- 
ing there should be smeared out. This makes extremely dif- 
ficult to constrain the innermost radius of the accretion disk 
from the iron-K emission profile. 



Ueda et al. 
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We find that a thermal Comptonization is the most likely 
origin of the continuum emission of GRS 1915+105 in "soft 
state". The standard, MCD plus a power law description for 
the continuum does not give acceptable fit. Instead of the 
MCD model, the p-fr ee disk model, expected from a one- 
dimensional slim disk (lAbramowicz et alii 19881) . gives a sta- 
tistically acceptable fit to the spectrum in Epoch JT . We con- 
sider it physically unlikely, however, because the inferred in- 
nermost radh«j^Joosmall ( 0.3r g ) than a theoretical pre- 
diction ( Watatai et al. 200(J. Also, we cannot consistently ex- 
plain the spectra of two epochs with the same model. 

The spectral model of Comptonization yields an optical 
depth « 5 and a low electron temperature («4 keV) for ei- 
ther of the two brightness levels that were sampled with the 
PCA. The seed photon spectrum is consistent with the emis- 
sion from a standard disk (MCD). No direct disk component 
is required from the fit. In this interpretation, most of the 
radiation energy is produced in the Comptonization corona, 
which completely covers the inner part o f the disk. Our re - 
sults basically support the conclusion by iDone et alJ (j.2004), 
who analyzed 7 representative PCA spectra of GRS 19 15+105 
in various states and found that thermal Comptonization is a 
dominant spectral component in most (5 out of 7) cases. In 
their color-color diagram (see their Figure 4), our spectra are 
located at the positions of (1.23,0.45) and (1.45,0.52), which 
are closest to obs . 7 and 6 for Epochs F and IF, respectively. 
IDone et ail (120041) favor a MCD-dominant model for both obs. 
7 and 6 rather than a purely Comptonized spectrum, however. 
The reason for this discrepancy is unclear 10 . 

In the Comptonization model, we can infer the innermost 
radius (r m ) of the accretion disk injecting the seed photons 
to the corona, by the conservation of photon number in the 
Comptonizat ion process. For this purpose , we adopt the for- 
mula (Al) of iKubota & Makishimal (|2004h with the left-hand 
term increased by a factor of two, 

/) r}„ cos(/) 7L , 19 

2^2cos(0 = 0.0165[^ I ^]( T ^) 3 photonss"' cm"", 

where F t £ c is 0.01-100 keV photon flux from the Comp- 
tonized component. Here we have assumed that half of pho- 
tons in the corona is injected again into the accretion disk due 
to the large optical depth (« 5). The results are given in Ta- 
ble 4. For isotropic emission of the Comptonized component, 
we find n n sal 10 km and «80 km for Epoch F and IF, respec- 
tively (instead, if the optically thick corona has a thin disk-like 
geometry, the 2cos(/) factor in the left term must be ignored 
and hence r m is increased by about 20%). The estimated r m 
values correspond to ss 5r g (Epoch F) and rj 4r g (Epoch IF). 
With correction for the color-to-effec tive temperature and that 
for the inner boundary condition (see Kubot aet al.ll 1998b . the 
physical radius would be larger than these values by a factor 
of 1.2. Thus, based on the innermost temperature and pho- 
ton numbers, the optically-thick disk may extend down only 
to £j 4r g in this state, not to the innermost stable circular or- 
bit (ISCO) of 1.23r g for a maximally rotating Kerr hole. If 
the radius observed in Epoch IF indeed corresponds to the 
ISCO, then our result indicates that the black hole spin of 
G RS 1915+105 is only m oderate, supporting the conclusion 
bv lMiddleton et al l (l2006l) . 

10 Done et al. (2004) introduce a strong reflection component (f2 re o/27r > 
2) in obs. 7 to explain the prominent edge feature around 9 keV as seen in 
our Epoch F spectrum. We infer that it may be significantly affected by deep 
absorption K-edges from highly ionized Fe ions in the disk wind, which are 
not included in their spectral model. 



Of all the classes of light curves for GRS 1915+105, the 
<f> class signifies conditions when the source is most soft and 
steady. In most BHBs, these conditions are associated with 
the high/soft state, in which a MCD component dominates 
the X -ray spectrum while there are lo w levels of Comptoniza- 
tion (iRemillard & McClintockl 120061) . normally observed at 
lower fractions of the Eddington rate than in our case (0.3- 
1 .2 L/Lsm depending on the geometry of the emitting region). 
Our analyses show that this is not the case for our <f> state ob- 
servations of GRS 1915+105, where Comptonization domi- 
nates. Furthermore, the 67 Hz X-ray QPO can be considered 
as additional signs of disturbances that lead away from the 
simple MCD that we envision for the high/soft state. Such 
disturbances are chronic prop e rties o f GRS 1915+105 in its 
soft states. iMcClintock et al.l (120061) were able to find only 
20 high/soft state examples in an archive of ~ 1000 RXTE 
pointings that included ^150 observations of soft and steady 
conditions (classes <fi, S, and some 7 types). 

Such strongly Comptonized spectra are characteristics of 
the very high state (VHS) o f other BHBs observed at lu mi- 
nosities of ~ 0.5L/LEdd (e.g jKubota & Makishimal2004l) . al- 
though the electron temperature of the Comptonizing corona 
in the VHS is much higher (> 10 keV) than our case («4 
keV). This suggests that the "soft state" of GRS 1915+105 
corresponds to the VHS of canonical BHBs, but the electrons 
in the corona must be more efficiently cooled by Comptoniza- 
tion. It could be explained because the disk extends to small 
radii in the "soft state" of GRS 1915+105, hence providing 
many soft photons to the corona, while the di sk of other BHBs 
in the VHS ma y be truncated at larg er radii ( Ku bota & Pond 
l2004t see also lDone & Ku bota 2006). To confirm this picture, 
more systematic studies of GRS 1915+105 in similar states at 
different luminosities are important. 

There is a puzzling combination of a fast QPO and spectral 
indications of strong Co mptonization, and this has been s een 
in other sources as well (Remillard & McClintock 2006). If 
the oscillations originate in the disk, then the Comptonizing 
corona must be sufficiently compact to avoid suppression of 
the 67 Hz oscillations due to scattering effects; for r e =5 the 
size must be much smaller than ~ 10r„. Alternatively, the 
formation of the corona may involve asymmetries or waves 
related to energy injection mechanisms. There is insufficient 
data to determine if all of the soft-state properties of GRS 
1915+105 that deviate from the high/soft state appear to be 
correlated. However, such knowledge would assist efforts to 
understand the cause of Comptonization and the mechanisms 
and interpretive values of the QPO frequencies. 
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